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Abstract 

We have statistically characterized the self-assembly of multi-layer polystyrene colloidal crys- 
tals, using the technique of vertical deposition, with parameters chosen to produce thick layers of 
self-assembled crystals in one deposition step. In addition, using a lithographically directed self- 
assembly method, we have shown that the size of multi-layer, continuous crack-free domains 
in lithographically defined areas can be many times larger than in the surrounding areas. In a 
single deposition step, we have produced continuous colloidal crystal films of 260 nm diameter 
polystyrene spheres approximately 30-40 layers thick, with a controllable lateral size of 80 - 
100 |Ltm without lithography, and as high as 250 \im with the lithographic template. 

Keywords: self-assembly, colloidal crystal, phononics, polystyrene spheres 



Email address: tero . isotalo® j yu. fi (Tero J. Isotalo) 
Preprint submitted to Journal of Colloid and Interface Science 



February 22, 2013 



1. Introduction 

Colloidal crystallization of mono-disperse spherical sub-micron particles has been used suc- 
cessfully in the fabrication of three-dimensional (3D) photonic crystals (TJ for a wide variety 
of applications [3]. These artificial opals, as they are commonly known, can also exhibit an 
acoustic or a phononic band gap, whose frequency depends on the particle diameter. While the 
properties of 3D phononic crystals have recently been studied by various groups in the millimeter 
and micrometer scale at acoustic frequencies |4, 5 , 6], much less work has been performed on 3D 
hypersonic phononic crystals with 100 nm scale dimensions, which produce phononic band gaps 
in the GHz range |7|. This frequency range is even relevant for thermal properties, as the dom- 
inant thermal phonons at low temperatures below IK have frequencies in that same range 181 . 
One possibility for GHz 3D phononic crystal fabrication is therefore the use of polystyrene (PS) 
nano-spheres of ~ 100 nm size range, as they are commercially available as colloidal suspen- 
sions. As demonstrated before (H [9j [TUJ [TT], these particles can be self-assembled into crystals 
from suspensions, allowing a convenient bottom-up assembly of a GHz phononic crystal with 
possible applications in thermal management or phononic waveguiding IT2l . 

Common methods for self-assembly of colloidal particles include gravitational sedimenta- 
tion, capillary assembly and vertical deposition 021. Typically, closed-packed lattices (face- 
centered cubic, hexagonal closed packed, or a random combination of the two) are produced in 
each of these methods, with face-centered cubic being more common. Here, we use the verti- 
cal deposition method for producing large area, thick 3D closed packed lattices using colloidal 
suspensions of PS nanospheres. This technique in its simplest form has been documented in the 
literature d [9l [TUl [TTJ quite extensively. Here, we concentrate on increasing the lateral single 
crystal domain size for the crystals and characterize the domain size results statistically, as a 
function of dipping speed and nanosphere concentration. We also combine the verical deposi- 
tion with lithographically patterned substrates, not to control the crystal structure itself as has 
been suggested before (HlttSHSlE), but to increase the size of single crystalline domains and 
to reduce cracking. In contrast, most other lithography directed methods (T5j [T6J [XT) are typi- 
cally mono-layer techniques. These have the drawbacks of multiple process steps and lengthy 
fabrication time, if thicker crystals are required. 

2. Material and methods 

Samples were coated with a solution of PS nanospheres using vertical deposition [ 1 ] and 
the resulting crystals were characterized by scanning electron microscopy (SEM). All substrates 
used in the experiments were either oxidized or nitridized silicon chips with dimensions of ap- 
proximately 8x16 mm, and coated with a layer of titanium oxide (TiO x ) approximately 10 nm 
thick on the surface. A titanium coating was first applied by electron-beam evaporation in an 
UHV chamber with a pressure 10" 8 mbar, and the resulting film was then oxidized in a vacuum 
chamber at 100 mbar of oxygen pressure for 2 minutes. The TiO x provided a hydrophilic sur- 
face for improved wetting by the nanosphere solution fT8l . A schematic of the dipping setup 
is shown in Figure [T] Samples without lithographically patterned surfaces were submerged into 
and withdrawn from solutions with concentrations of 0.02%, 0.2%, 2%, 5% and 10%. The orig- 
inal solution of 10% concentration, purchased from Duke Scientific, was diluted with de-ionized 
water to produce the lower concentrations. Each concentration was tested at withdrawal speeds 
from 0.01 mm/min up to 0.04 mm/min. It was determined that the speed of withdrawal from 
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Figure I: Schematics of the dipping setups are shown. TiO x coated substrates were dipped 
vertically and also at a 45° angle. 



the solution was the dominant factor in self-assembly and thus, a higher speed was used for sub- 
mersion in order to reduce total dipping time. In addition to vertical dipping, some experiments 
were performed at a 45° angle for select combinations of dipping speed and solution concentra- 
tion. Lithographically assisted self-assembly was investigated with patterned surfaces made by 




Figure 2: Lithography process steps are shown, (a) Pattern is exposed in e-beam resist (PMMA), 
(b) resist pattern is developed, (c) SiN x is etched with CF 4 plasma and then (d) the Si substrate 
is etched anisotropically with KOH. (e) Second lithography is aligned with the previous pattern, 
(f) pattern is developed, (g) UHV evaporation of Ti and oxidation followed by (h) final lift-off. 

a two-step electron-beam lithography technique, shown scematically in Fig. 2. Samples were 
fabricated by first etching troughs into silicon substrates, using the overlay SiN x layer (750 nm) 
as an etch mask. The nitride mask was defined by etching it reactively in CF4 or CHF3 plasma 
(for CHF 3 , we typically had to etch for 20 min with parameters min 100 W, 55 mTorr). The 
troughs with sidewalls angled at 54.74° were produced by wet chemical etching in potassium 
hydroxide (KOH) at 90 °C. The depth of the troughs were varied from a few \im up to nearly 
100 |Ltm. These troughs were then selectively coated with a hydrophilic TiO x layer by a second 
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e-beam lithography step with litf-off aligned with the previous pattern, which produced TiO x 
coating only at the bottoms of the troughs. 



3. Results 



The fabricated crystals typically consisted of a set of domains, separated by cracks produced 
by the drying process. Domain size data were collected by taking several SEM images from 
each sample, and using the scale bar as a calibration length. A representative measurement 
image is shown in Figure [3jleft). All measurements take into account the direction of dipping, 
where measurements taken along the dipping direction were defined as vertical lengths, while 
those which are perpendicular were defined as horizontal. Also shown in Fig. 3 (right), is a 
side-view close up of a typical thick multi-layer colloidal crystal structure. Large data sets were 
collected from each sample in order to improve statistics and reduce the inherent roughness of 
the measurement method. From the scatter plots of vertical vs. horizontal length, we always saw 
a distribution of domain sizes for all speeds and solution concentrations. 




Figure 3: An example of a SEM image used for the domain length measurements (left,) with a 
side-view close up of the multi-layer colloidal crystal structure (right). Values shown for mea- 
sured lengths are in |Ltm. 



The plain, unpatterned samples dipped into solutions of 0.02% and 0.2% concentration did 
not produce any crystalline order even at the lowest dipping speeds and thus, no data is shown 
for these concentrations. For the higher concentrations crystalline domains always formed, and 
an example of the scatter plot of domain size measurement results for 2%, 5% and 10% solutions 



dipped at 0.02 mm/min is shown in Figure 4a From these, a clear trend can be seen with 



increasing concentration. While there is a trend toward larger average domains with increasing 
concentration, the size distribution also increases. This variability in domain size becomes quite 
large, as seen in the data for 10% concentration, where the vertical domain size varies from a few 
tens of |Ltm up to approximately 150 |Ltm. Similar trends are seen in the dipping speed comparison 
plots, with an example shown in Figure 4b for the 10% solution, at four dipping speeds. At 0.01 
mm/min, we see vertical lengths ranging from several tens of m up to approximately 350 |Ltm. 
Again, while the domain size grows with slower dipping speed, so does the size distribution. 
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Dipping Speed 0.02 mm/min 1 % Concentration 




Horizontal Length (urn) Horizontal Length (um) 



Figure 4: (a) A comparison of domain sizes is shown for concentrations of 10%, 5% and 2%, 
with all concentrations dipped at 0.02 mm/min, or (b) for dipping speeds between 0.01 mm/min 
and 0.04 mm/min. The maximum domain size as well as the distribution increase with decreasing 
speed and increasing concentration. 
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Figure 5: (a) Average vertical length of domains is plotted against dipping speed for vertical and 
angled dipping. Clear trends can be seen for dipping speed, concentration and angle, (b) Ratios 
of vertical to horizontal length of domains are plotted for vertical dipping. Vertical domain 
elongation and general spreading of domain sizes can be seen. 
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Figure |5a| shows the average vertical domain size for vertical and 45° angled dipping, using 
three concentrations at four dipping speeds. Angled dipping was found to produce smaller do- 
mains on average. While the exact cause of this is unclear, it may be that changes in the shape 
of the meniscus can alter the evaporation rate at the self-assembly region. This could lead to 
shorter effective self-assembly times in the angled dipping cases. Though trends are clear, the 
wide distributions in domain sizes result in some scatter. For example, at 0.02 and 0.03 mm/min 
vertical dipping speeds, we see that the average domain size does not correlate perfectly with 
concentration. In the angled dipping, the 5% concentration shows negligible change in average 



domain size with increasing speed. While Figure 5a shows an apparent asymptotic approach to 
some minimum domain size, such a minimum is not supported by the data from the lower con- 
centration solutions, where no domains were found. Additionally, an indication of directionality 
in domain growth can be seen from ratios of the vertical to horizontal length. This elongation ef- 



fect can be seen from Figure 5b which shows these ratios for three concentrations at four dipping 
speeds. (The positions of data for 2% and 10% concentrations are shifted in the plot for visual 
clarity). The average ratio, shown in red for each data set, indicates an increased elongation 
along dipping direction with increasing dipping speed. These data also show that the distribution 
of ratios widens with increasing dipping speed. 
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Figure 6: Histograms of vertical domain lengths are shown for (a) 2%, (b) 5% and (c) 10% 
concentrations at a dipping speed of 0.01 mm/min. All three concentrations are characterized by 
a log-normal size distribution. 



The domain size distributions for all dipping speeds and concentrations were also investi- 
gated, with example histograms of vertical domain lengths shown in Figure [6] for 2%, 5% and 
10% concentrations at 0.01 mm/min dipping speed. We found that all concentrations were ap- 
proximately characterized by log-normal size distributions shown in red, especially at the lowest 
dipping speeds. Only at the highest dipping speeds do the distributions start to be less skewed 
and approach the normal distribution. Thus, most of the data are not well described by arithmetic 
mean and standard deviation. Compiling, instead, the geometric means x g = (H n x n ) l/n and the 

geometric standard deviations cr g - exp ^ ^2^(ln 2 (x„/Xg)/ftj of these distributions, we found 

that all cr g values fell in the range of 1.32 and 1.71, with the exception of two anomalously wide 
distributions (at 2% concentration and 0.02 mm/min dipping speed, the horizontal and vertical 
domain size distributions had deviations cr g of 2.025 and 2.017, respectively). The normalized 
geometric standard deviations cr g /x g are very well approximated by a l/x g trend for nearly all 
data points, as shown in Figure [7] except for the two anomalous cr g /x g values seen (circled in 
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Figure 7: Normalized geometric standard deviation cr g /x g data for all dipping speeds and con- 
centrations vs. geometric mean. Two anomalously wide distributions (circled in red) can be seen 
to fall farther from the 1 /x fit (dashed curve). 



red) in Figure [7] slightly farther from the fit curve. These two sets of data were likely the result 
of insufficient data sampling or poor substrate quality. 

The effects of etched lithographic features on self-assembly were studied by comparing re- 
sults of samples dipped at various speeds with troughs of varying size and depth. For troughs 
of 10 |Ltm depth, there was hardly any change in domain sizes compared to areas outside the 
trough, while troughs of 30 \im depth did not appear to fill completely. All troughs of 20 \im 
depth showed more consistent order in general than unpatterned areas. The lateral shape of the 
trough was always a square, but with a side length varying between 50 - 100 \xm in steps of 10 
|Ltm, and for larger troughs from 150 to 300 \im in steps of 50 \im. Troughs smaller than 70 |Ltm 
showed little to no difference in domain size when compared to the surrounding area. In troughs 
larger than around 200 |Lim, there were always more than one domain. Figure 8a shows clearly 
the increase of domain size within troughs compared to the unpatterned areas. Domain size out- 
side the trough is approximately 30 - 70 |Lim while inside, the crack-free domains can be many 
times larger. The consistency of the continuous crystals in troughs indicates that this method can 
be used to locally suppress the domain size distribution completely. Figure 8b shows a closer 
view of the largest trough and its surroundings. While crack- free domains up to approximately 
200 \im size are reproducible, these regions are not necessarily single crystalline. All samples 
produced inside the troughs so far have polycrystalline order of crystallites with differing crystal 
orientations, and the size of these single crystal domains -10 |Lim is consistent throughout all 
the samples. 
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(a) (b) 




Figure 8: (a) SEM image of a sample with a range of etch trough sizes, after dipping in 10 % 
solution with speed 0.02 mm/min. Areas outside the troughs have domains ranging from 30 - 
70 |Ltm while the etched troughs show significantly larger domains, (b) Close up view of a 250 
|Ltm by 250 |Ltm trough, seen with a continuous crystalline region several times larger than the 
surrounding domains. 

4. Conclusions 

Vertical deposition on TiO x coated substrates has been used to produce self-assembled col- 
loidal crystals of polystyrene nano-spheres. For unaltered surfaces, using a solution of 10% PS 
spheres and a dipping speed of 0.01 mm/min, we achieved an average vertical domain size of 
approximately 130 |Ltm with some domains up to 350 \xm. These unaltered substrates produce, 
however, large domain size variations, characterized by log-normal distributions, tending only 
slightly toward normal distributions at higher dipping speeds. The relative variation in domain 
size increases with both dipping speed and concentration. Nevertheless, the variability even at 
the lowest speeds is significant and poses a serious challenge for the fabrication of photonic or 
phononic crystals on largest scales. 

On the other hand, lithographically assisted self-assembly was used to consistently produce 
continuous regions of multi-layer colloidal crystal as large as 200 \xm, effectively suppressing 
the size distributions found in unaltered substrates. While these regions are free of cracks, they 
are still not fully single crystalline, and the crystalline quality of these lithographically assisted 
multi-layer structures was comparable to that of the plain self-assembly samples. Improving 
the crystalline quality of these regions remains a significant challenge for rapid production of 
large-scale structures. Still, the method investigated in this paper is notably faster than most 
methods for large-scale colloidal crystal fabrication and thus, may be of interest to the photonic 
and phononic crystal communities. 



8 



References 



[1] P. Jiang, J. F. Bertone, K. S. Hwang, V. L. Colvin, Chemistry of Materials 11 (1999) 2132-2140. doi |10 . 102l7| 
cm990080+ 

[2] Y. A. Vlasov, X.-Z. Bo, J. C. Sturm, D. J. Norris, Nature 414 (2001) 289 - 293. doi[l] 

[3] J. D. Joa nnopoulos, P. R. Villeneu ve, S. Fan, Nature 386 (1997) 143-149. URL: ]http : //www . nature . com/| 

Idoif inder/10. 10 38/386143a0l 

[4] Z. Liu, X. Zhang, Y. Mao, Y Y Zhu, Z. Yang, C. T. Chan, P. Sheng, Science 289 (2000) 1734-1736. doi |10 . 1126/] 

[science . 289 . 548 5. 1734] 

[5] S. Yang, J. H. Page, Z. Liu, M. L. Cowan, C. T. Chan, P. Sheng Ph ys. Rev. Le tt. 88 (2002) 10 4301. URL: |http : | 

|//link . aps . org/doi/10 . 1 10 3/PhysRev Lett . 88 . 1 04301 dor 10 . 1103/PhysRevLett . 88 . 104301 
[6] S. Yang, J. Page, Z. Liu, M. Cowan, C. Chan, P. Sheng, Physi cal Review Letters 93 (2004) 2-5. URL: |http H 

//link . aps . org/doi/10 . 1 103/PhysRevLett . 93 . 024301 
[7] W. Cheng, J. Wang, U. Jonas, G. Fytas, N. Stefanou, Nature Materials 5 (2006) 830 - 836. doi 1 
[8] N. Zen, T. A. Puurtinen, T. J. Isotalo, I. J. Maasilta, submitted (2013). doi[l] 
[9] Z.-Z. Gu, A. Fujishima, O. Sato, Chemistry of Materials 14 (2013) 760-765. doi[l] 
[10] C.-A. Fustin, G. Glasser, H. W. Spiess, U. Jonas, Langmuir 20 (2004) 9114-9123. doifT] 

[11] R. G. Shim min, A. J. DiMauro, P . V. Braun, Langmuir 22 ( 2006) 6507-6513. URL: |http : //pubs . acs . org/| 

|doi/abs/10. 1021/la060278q| doi |10 . 1021/la 060278q pMID: 16830991. 
[12] Y. Pennec, J. O. Vasseur, B. Djafari-Rouhani, L. Dobrzynski, P. A. Deymier, Surface Science Reports 65 (2010) 

229-291. URL: http://w ww.sciencedirect. com/scie nce/article/pii/S 0167572910000555] doi |10.| 

|1016/ j . surf rep . 2010 . 08 . OO^T 
[13] Y Xia, B. Gates, Y Yin, Y Lu, Advanced Materials 12 (2000) 693-713. doi 1 
[14] A. van Blaaderen, R. Ruel, P. Wiltzius, Nature 385 (1997) 321-324. doi[l] 

[15] Y. Xia, Y. Yin, Y. Lu, J . McLellan, Advanced Functional Mate rials 1 3 (2003) 907-918. URL: |http://dx.doi.| 

|org/10.1002/adfm. 200300002 doi | 10 . 1002/adfm. 200300002 

[16] N. N. Khanh, K. B. Yoon, Journal of the American Ch emical Society 131 (20 09) 14228-14230. URL: jhttp :"| 

|//pubs . acs . org/doi/abs/10 . 1021/ j a905534k] doi |10 . 1021 /ja905534k pMID: 19775130. 
[17] J. Zhang, Y. Li, X. Zhang , B. Yang, Advanced Materials 22 (2010) 4249-4269. URL: |http : //dx . doi . org/10 . | 

|1002/adma. 2010007551 doi |10 . 1002/adma. 201000755| 

[18] S. Ka ruppucham y, J. M. Jeong, Materials Chemistry and Physics 93 ( 200 5) 251 - 254. URL: |http: //www 7 

|scie ncedirect . com/science/article/pii/S0254058405002579 doi 10 . 1016/ j .matchemphys . 2005 . 

lo4.015l 



9 



